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Abstract 



.^^ ' The dual-tree complex wavelet transform (DT-CWT) is an enhancement of the con- 

ventional discrete wavelet transform (DWT) due to a higher degree of shift-invariance 
and a greater directional selectivity, finding its applications in signal and image process- 
r^ ' ing. This paper presents a quantitative proof of the superiority of the DT-CWT over 

the DWT in case of modulated wavelets. 
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>. ! 1 Introduction 

(N 

en I Wavelet transforms provide a convenient technique to perform a multiresolution analysis 

O^ . of finite-energy signals. The most popular instance of a wavelet transform is the critically 

^ . ' sampled discrete wavelet transform (DWT), which is an invertible transform that permits 

'nI" . sparse signal decompositions at a low computational cost [T^ . 

^—^ I The DWT has been successfully employed in many applications, including image compres- 

sion [20], noise reduction [8j and speech recognition |9J. However, in the area of statistical 
signal processing, the DWT has proven to be less effective [3 [T3]. This is mainly due to 
the high translation sensitivity of the DWT: small shifts in the input signal may completely 
^^ ' change the wavelet coefficient pattern. As a consequence, algorithms based on the DWT need 

^ . to recognize and understand a wide variety of different wavelet patterns. 

One way to address the shift- variance problem is to relax the critical sampling criterion 
of the DWT. In [TS], an overcomplete version of the DWT is proposed, which is most easily 
implemented by the "a trous" algorithm. A generalization of this algorithm is described in [3] . 
Note that this approach is computationally intensive and produces highly redundant output 
information, which limits its applicability. Nevertheless, since the output of the "a trous" 
algorithm can be computed directly from the critically sampled DWT, it is readily applied in 
DWT-based image and video coding systems \T, T' . 

In [181I19| . Simoncelli et al. introduce the steerable pyramid, an alternative decomposition 
method based on Laplacian pyramids and steerable filters that achieves approximate shift- 
invariance. Furthermore, the steerable pyramid also gives a better directional selectivity 
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when analyzing two-dimensional signals, which simplifies the extraction of geometric features 
in images. 

Another way to improve the shiftability of the DWT is by simultaneously employing two 
real DWT channels that form an approximate Hilbert Transform pair. By combining the 
corresponding coefficients of the first and second DWT into complex- valued coefficients, we 
obtain a new transform, which is called the dual-tree complex wavelet transform (DT-CWT). 
Compared to the steerable pyramid, the DT-CWT provides a better directional selectivity 
while having a lower redundancy factor of 2'^ for d-dimensional signals. A more elaborate 
discussion on the design and use of the DT-CWT can be found in [17] . 

The near shift-invariance property of the DT-CWT has been extensively studied over 
the last decade [51 [TUl [HI [131 [HI [23 ■ Recently, Chaudhury and Unser [5j deduced an 
amplitude-phase representation for dual-tree complex wavelet transforms that involve m,od- 
ulated wavelets, linking the multiresolution framework of the wavelet components to the fre- 
quency decomposition through Fourier Analysis. This representation provided new insights 
into the improved shiftability of the DT-CWT. 

In this contribution, we build on their findings by introducing a more formal description 
of the DWT translation sensitivity, which will allow us to better explain the superiority of 
the DT-CWT. We finish with a study on the decaying rate of the DT-CWT shift error when 
the translation parameter tends to zero in case of orthonormal wavelet systems. 

2 Preliminaries 

We now introduce some definitions and notation needed to state our results in the following 
sections. 

Given two signals / and g in L^(]R), we define their inner product by 



(/,.9) = / f{x)g{x)dx, 
where the bar indicates complex conjugation. The Fourier transform of / is given by 

/(O = / f{x)e~^^^dx 

Jr 

whereas the Hilbert transform H (HT) is characterized by the relation 

The Hilbert transform is orthogonal to the signal, commutes with translations and positive 
dilatations, and Ti^^ = —T-L. The translation-dilatation operator Sj ^ on -0 e L^(IR) is defined 

by 

Let {tlij,k\j,kez and {ip'j k}j,kez be two real-valued bi-orthogonal wavelet systems that 
form a Hilbert transform pair, i.e. ip' — 'Hip. We define the wavelet coefficients of / with 
respect to these wavelet systems by 

ajW = (/>j\fc) and bj[k] ^ {f,^}'^^) 

for every j, k E Z. These equations yield the following two different wavelet identities: 

/= XI «iW^J\'= and f ^ Yl ^iW^i.fc' 
j,ke'i j,kei 

where ij^j^k and ■0^ ^ represent the dual wavelets of il^j.k and i/)' j. respectively. 
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for every j, fc G Z. 

Recall that for dyadic wavelet transforms, the level j coefficients of a shifted signal /(• + s) 
with s = 2~^m, m S Z, can be easily predicted from the coefficients of the reference signal. 
In fact, 

c,^[fc] = (/(. + 2--'m), *,,,.) 
= (/,%fe(.-2--'m)) 

= Cj [fc + to] . 

This well-known property can be adapted for arbitrary shifts s by decomposing s into a dyadic 
number 2~^m and some remainder h with \h\ < 2^-': 

s = 2^^m + h. 



Then 



c. 






= C^[fc + TO]. 



The adjusted shift error 



\cj[k + to] - cj[fc]| = |cj[fc + to] - c^[fc + to]| 

is in general much smaller than the original shift error \cj[k] — c^[fc]|. In order to further 
reduce the shift error \cj[k + to.] — c^ [fc + to,] | , we will perform a phase change of Cj[k + in] 
over an angle (ph that partially compensates for the small shift h {\h\ < 2~^), so that 

cj[fc"] = c'^[k + to] « e"'"^Cj[k + to]. 

As suggested in [S], we make the assumption that the involved wavelet ip is modulated. 
That is, 

ip{x) = w{x) cos(ajoa; + ■^o) 

for Wo, Co > where the localization window w is bandlimitcd to [—0,17] for some D, < ujq. 
Examples of modulated wavelets are the Shannon and Gabor wavelets. As the orthonormal 
spline, resp. S-spline, wavelets resemble the Shannon, resp. Gabor, wavelet, they can be 
seen as a kind of modulated wavelets. Using the Bedrosian identity (see [5]), one can show 
that 

ip'{x) = w{x) sin(ajoa; + Co)- 



In this way, we obtain the identity 

In order to examine the near shift-invariance of the DT-CWT based on these modulated 
wavelets, we thus need to minimize the error 



\e^^-c,[k]-c'j[k]\ 

for some well-chosen angle (j^h S ["""i i'[- This phase- compensated shift error will be compared 
to the "real" shift errors \aj[k] — a^[/i;]| and \bj[k] — b'^[k]\ in Section [3l The attained results 
(summarized in Theorem l3.1|) suggest to take 4'h = 2^ujQh, which accords with the conclusions 
drawn by Chaudhury in 6, p. 127-128]. 

/,8 



We have empirically verified that c^[fc] « e'^^"°''cj[/c] for small shifts h using the DT-CWT 
software provided by the authors of [4]. In our experiments, we put ip = ''Pt 5 ^nd '0' = f/if, 
the fractional B-spline wavelets of degree a = 8 and with shift parameters r = 4.5, t = 5 
respectively [21j . These wavelets are known to be approximately modulated. In fact, one can 
observe that 

i{}{x) « w{x) cos(5.3x -I- 5.2), 



where w ~ ^/ip^ + ^'^. It is proven in [3] that ip and ijj' form a Hilbert Transform pair; 
therefore, they determine a DT-CWT. 

Our test signal consists of 512 uniform samples from a block function / : [0, 1] — !> M (see 
Figure [Ij, which is shifted one place to the left. The resulting signal corresponds to the 
function f^ = /(•-!- 1/512). 

Figure[2]compares the phase-compensated shift error |e'^ "°''cj[fc]— c^[fc]| with the optimal 
shift error 

The graphs on the right are zoomed out so that the shift errors of the real and imaginary 
wavelet components can be included. The following naming conventions are used for the 
considered shift-errors: 

complex-optimal = |cj[A:]| — |c[fc]| 
complex-phasecomp = cxp(i * 2-' * 5.3/512)cj[fc] — c [k] 
real = 2|Re(cj-[fc]) - Re{c'j[k])\ 
imag = 2|Im(cj[/c]) - Im(Cj^[A:])| 

The test results show that the phase-compensated shift error is significantly smaller than 
the shift error of the real and imaginary wavelet components. Moreover, there is a very good 
fit between the phase-compensated and the optimal shift error. Therefore, we can conclude 
that the derived phase-compensation for the DT-CWT coefficients is near-optimal under small 
translations. 



Test Signal 



B-spline HT pair with a=8 and t=(4.5,5) 




1 128 256 384 512 



Figure 1: Graphical representation of the test signal and the selected fractional i?-spline 
wavelet pair. 
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Figure 2: Graphs showing the shift-errors of the DT-CWT and the real and imaginary wavelet 
components. The plots on the left compare the phase-compensated error with the optimal 
shift-error. The plots on the right are zoomed out so that they can show the shift errors of 
the real and imaginary wavelet components. 



3 On the phase-compensated shift error of the Dual- Tree 
Complex Wavelet Transform 

Let / be a real-valued function in L^(R) with DT-CWT coefScients Cj[k],j,k G Z, based on 
two modulated wavelets ip and f/'' — Tiip. Denote the translates of / over some real number 
hhy f^ = /(• + h). The DT-CWT coefficients of /^ are given by 

As deduced in Section [21 to study the shift error for the level j coefficients, it is enough to 
look at \h\ <2~^. 

In this section, we work towards a proof of the following theorem. 

Theorem 3.1 (Phase Compensation for the DT-CWT shift error). Let Cj[k], j,k G Z, be 
the coefficients of a real-valued function f in L'^(R) with respect to a DT-CWT decomposition 
for which the involved wavelet -0 is modulated. When h is small, we have the approximate 
identities 

|e"'-o^c,[fc]-c^'[fc]| \e^^'-'-^c,[k]-c^[k]\ 

~' and ~ '^ 



\a,[k]~a>^[k]\ \h,[k]-h>^[k]\ 

proving that the phase-compensated error is negligible in relation to the shift errors of the real 
and imaginary wavelet components. 

We express the translation sensitivity of the DWT coefficients Oj [k] and bj [k] by postu- 
lating that 

\aj[k]-a^Ak]\ \b j[k] - b)[k]\ 

mf — -f = Ba and nif ., , ..■',. = Bt, (1) 

\h\<2-3 \ha-j[k\\ |/i|<2-j |/i6j[fc]| 

for some values Ba and Bi, significantly larger than zero. 

In the next proposition, we introduce the ratio Rfn which relates the phase-compensated 
shift error |e"^''Cj[/c] — c'NA:]| to the shift errors of the real and imaginary wavelet components 

K[fc]-a^^[fc]|and|6,[fc]-6,^[fc]|. 

Proposition 3.2. Let (f)h G [— tt, 7r[ for every h E M. with \h\ < 2^-'. Define 

_ e^1'^a,[k]-a'}[k] 
'' e'^>^bj[k] - b'^ikY ^ ' 

Suppose that the constraint in (pp holds for Ba and B}, significantly larger than zero. Then 



3 I 

and 



\e^*^c,[k]-c'^[k] \ ^^^ ^^^^\l-^/R^\ 
\a,[k]-am\ 



<il + <P/Ba y^ r''"' (3) 



''^c,\k]-^m\^,^ ,...jfi.-^i 



\b,[k]-b';[k]\ 

where $ = sup|,j|<2-j |(e*"*'' - 1)/^|- 

Note that $ is finite if and only j/limsup^^Q \(j)h/h\ < cxd. 



,u.i. <a + '^/Bb)^^^^, (4) 



Proof. We prove ([3]); the proof of (U) is similar. Observe that 



= ^\l-i/Rh\\e'*'^a,[k]-a'j[k]\. 
Dividing \e"'"^Cj[k] ~ c'}[k]\ by \e^'t"^aj[k] ~ a!l[k]\ results in 



-J'*"^Cj[k] -c^;[k] 



On the other hand, we have 



which implies that 



|e''^''-l|KWI 
\a,[k]~am\ 



\l-t/Rh\ 



< '^/Ba, 



\e^^^aM - «,^[fc]l ^ h#] - a'jM + le'*" - 1| |a,[/c]| 



Hence, 



\a,[k]-a)[k]\ - \a,[k]-a^^[k]\ 

< 1 + ^/Ba. 

\e^'t>-c,[k] - c'^m _ \e^*'^a,[k] - a^'[/c]| |1 - ^/Rf,\ 



,[k]-am\ 



\a,[k]-am\ 



<(i + $/B„: 



\l-i/Rh\ 



D 



This shows that the phase-compensated shift error |e"^''Cj[/c] — c'j[k]\ becomes smaller as the 
ratio Rh approaches to i. 

The perturbed coefhcients c^ [k] can be expressed more explicitly as 



-i£,o 



f{x)Ej,k[wix - 2^7i)e-*"''(^-2''')]da;. 



Observe that both the sinusoid and the localization window of the modulated wavelet con- 
tribute to the perturbation of c^ [k] . Their individual roles on the DT-CWT shiftability can 



be described using the variables 



Eh^ I f{x)E,^k z«(x)(e-'^"«^-e'^"«("-2^")) 



and 



Wh 



f{x)^j^k {w{x) ~w{x -2^h))e 



dx 



dx. 



Proposition 3.3. If the localization window w is differentiahle, then \Wh/Eh\ converges as 
/i — ^ 0. More precisely, we have the identity 



1 



h^o 2ujQ\cj[k\\ 



f{x)E,M^{x)e-^-o^]dx 



By the Cauchy-Schwarz inequality, we get that 

1 



1^-?oI^''/^'^i^2^;;^mII^IMI^ 



2 1 1 dx 1 1 2 ■ 



This expression reveals three important parameters that influence the asymptotic behavior 
of \Wh/Eh\ as h ^ 0: firstly, the frequency wq of the modulated wavelet; secondly, the L^- 
norm of ^ ; and thirdly, the significance of the DT-CWT coefficient Trap- • 

The first two parameters are close to zero for most modulated wavelets, because the overall 
frequency ujq is very high compared to the slowly varying localization window w, which is for 
example illustrated in Figure [TJ We can thus conclude that \Wh/Eh\ « for all significant 
coefficients Cj[k] when h is small. 

Clearly, 

\E,\^2\e^''-'''^-l\\cM- (5) 

By defining ah as the unique angle in [0,7r[ that satisfies 

|W^,| = 2|e-'^-l||c,[fc]|, (6) 

we arrive at 

In the next proposition, we show that the ratio Rh corresponding to the phase-compensation 
(ph = 2^ujQh + sign(ft,)aft is approximately equal to i when h is small and \Wh/Eh\ ~ 0. 

Proposition 3.4. Let (3h g] — tt, tt], such that f3h = 2^^^u!Qh modulo 2-k. Define 

(j)h = 2'ujoh + sign{l3h)ah. (8) 

Ifah<TT~ IPhl and \Wh/Eh\ < I, then 

^ e^'l>-a,[k]-a'}[k] ^ . l + Kf, 
'' e^^-h,[k]-h^^[k] 'l^Kh' ^' 

where K^ is a complex number that satisfies 

2\Wh/Eh\ 



Kh < 



|gi2i<^o/i_^l|_ \Wh/Eh\' 



Proof. We prove the proposition for /3/i > 0. Define A — e"^'*/ — /''. Since V-" = ^ + ^ and 
-0' — — i(^ — \1/) by construction, we obtain that 

e'^'^a,[k] - a^'[A;] = (A, *,-,.) + (A, *j- ,-) 
and 

e'^'^h,[k] - h';[k] - i(A, *,■ fe) - i(A,'^j- fc). 
This shows that Rh = i(l + Kh)/{1 - Kh), with 

_ (A,^,,,) 



A simple calculation gives 

(A,*,,,)=e^^"(/,*,,,)-(/,vl/^-,.(.-/i)) 



-iio 



/(a;)Sj-,fc e-''^'"'[e"'"^wix) - e''^'"«^w{x - Vh)] 



^%{Vu,ah-io) 



-Wh + {e"'"^ -e'^''^°^)cj[k]. 



Similarly, 



Hence, 



(A, *,,,.) = Wh + (e"^^ 



\Kh\< 



^-i2'uoh\ 



[k]. 



Wh + 2 eJi't>h-2^u:oh) _ 


-1 c,[k] 


Wh -2 e'K'^h+2^^oh) _ 


-1 c,[k] 



The substitution (j)h = 2^ujQh + ah gives 

\Wh\ +2|e''"" -l| \cj[k]\ 



\Kh\< 



\Wh\ - 2|e*(2^+''^«''+"'-) - ll \ci\k]\ 



Since 



|g«(2^+ uoh+a^) _ i|2 = 2 - 2 cos{Ph + a/i) 
>2-2cosil3h) 

^ l^t2^ + ''uJoh _ 2|2 



we have 



< 



1 



gi(2J + iwo/i+Qh) - 1 ~ gi23 + iwo/!. _ 1 

1 



\p^i2^ujoh ^\ \p^i23ujQh i ]^ 



This gives us, in combination with Equation 

\Wh\ 



gz(2J + ia;o/i+ah) — ll ~ |gi23LJo'i + 11 '' 



On the other hand. 



We arrive at the bound 



Wh 



giafe _ 1 



<2|c,[fc]| 



\Kh\ < 



4\cM\Wh/Eh\ 



2|cj[/c]| |e^2^"t"' + l| - 2\cj[k]\ \Wh/Eh\ 
2\Wh/Eh\ 



\e^"-«'' + l\-\Wh/Eh\' 



This finishes the proof. 
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Figure 3: Plot of the ratio Rh for h — 1/512 and j — 3. Observe that Rh differs more from i 
for fc = 4, 5. 



The phase-compensation (j)^ proposed in (|8]) may be hard to determine in practice. Instead, 
when h is small and \Wh/Eh\ « 0, we can put 0/1 — 2^ujoh. Indeed, Formula ^ implies that 



lim \Wh/Eh\ = -^ lim ^. 



As a consequence, 



(I)h/Vi0oh - 1 + signih)ah/2^LUoh « 1 + \Wh/E,,\ w 1. 



(10) 



Figure |3] plots the ratio Rh for a number of DT-CWT coefficients, using the same test con- 
figuration as described in Section [2j We see that the values of Rh are indeed close to i. Also 
note that the accuracy of the approximation depends on the magnitude of the DT-CWT 
coefficients. This phenomenon is to be expected, as the size of \Wh/Eh\ increases when the 
significance of the coefficient Cj[k] decreases (see Proposition [ 



Proof of Theorem \3.1l Let 4>h = 2^ujQh. As noted after Proposition 13.31 we may assume that 
\Wh/Eh\ ~ 0. Then (jlOp shows that this instance of (f>h is equivalent to the one specified 
in ([8]). Hence, Proposition 13.41 is applicable, resulting in Rh ~ i- Substitution into ([3|) and 



Q of Proposition [321 reveals that the ratios 



|e"'^-°''c,[fc]-c^^[fc]| |e'"'"°''c,W-cj-[fc]| 



\a,[k]-a^[k]\ 

approximately zero. This is exactly what we needed to prove. 



and 



|fc,[fc]-6H/c] 



are both 
D 



4 On the decaying rate of the phase-compensated shift 
error 

In this section, we focus on the dual-tree complex wavelet transform for which the wavelet 
systems {V'j.fcjj.fcez and {ipj k}j,kei. are both orthonormal and modulated. In order to estimate 
the prediction error c^[fc] ~ e*^ '^'^^Cj\k]^ we will introduce a new quantitative bound for the 



phase-compensated shift error |e* 



, [k]—c^ [fc] I . Moreover, this bound allows us to describe 



the decaying rate of the phase-compensated shift error as /i — ?► 0. 
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In [S] , the fractional Hilhert Transform (fHT) operator is introduced in order to deduce an 
amplitude-phase representation of the DT-CWT. The fHT corresponding to the real- valued 
shift r is defined as 

Hr = cos(7rT)I — sin(7rT)?^, 

where I is the identity operator. Note that for t = —1/2, we retrieve the original Hilbert 
transform operator. Moreover, 'Ht[cos{ujqx)] = cos{ujox + ttt). It is easy to show that the fHT 
is a unitary operator that commutes with translations and positive dilatations. In particular, 
if {ipj,k}j,kei, is a wavelet system, then {7^r['0j,fc]}, ..^^ ^^ ^^^^ ^ wavelet system. 

Theorem 4.1 (Amplitude-phase representation of the DT-CWT, [5 ). Let f be a function 
in i2(R) mth DT-CWT coefficients Cj[k] = |cj[fc]|e'"^W. Then 

f=Y^ |c,[fc]|S,.fc[u;C(^,[fc])], (11) 

j.kez 

where, for lo G [0, 27r[, C(ui) is given by 

C{lu){x) = cos(a;oa; + Co + ^)- 

Using this theorem, the authors of [S] provided new insights on the shiftability of the 
DT-CWT. Indeed, Formula ((TT|) gives an explicit interpretation of the phase parameter luj [k] 
as the phase-shift applied to the modulated sinusoid of the wavelet. More precisely, when / 
is shifted over h, we get an amplitude-phase representation of the form 

j.kei, 

Hence, the localization window w is kept fixed at scale j while the oscillation is now shifted 
over u!j'[k] to better fit the underlying signal singularities/transitions. 

We now extend upon their findings by employing the phase-compensated shift error to 
characterize the shift errors |e*^ ""''cj [k] — rf [fc] | , as stated in the next proposition. This result 

will be applied in Theorem 14.31 to estimate the prediction error c^[k] « e^^^'^°'^Cj[k]. 

Proposition 4.2. Let f be a real-valued function in L~ (M.) with DT-CWT coefficients Cj[k] ~ 

|Cj[fc]|e*'^^W. Consider a translate f' = fi-+h) with DT-CWT coefficients c'f[k] = | c^ [/c] | e*"^^" ^ . 
Then 



z'^''^'>''cj[k]-c'l[k]f 



je.hkeK 



where J, K '^ Z, 



ei 



Y, IcjM 2„fc [{w - w{- + Vh))CXu,[k] + 2^LUoh)] 
jeJ,keK 



and 



£2 



^ \cj[k]\ S,, fc [{w - w{- + 2^h))C{-uj,[k] - 2'LJoh)] 
jeJ,keK 
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Proof. Let us first prove the theorem for J = X = Z. Recall that Cj[k] = ^{aj[k] — ibj[k]), 
where Qj [k] and bj [k] are the coefficients of / corresponding to the real and imaginary wavelet 
components respectively. The same formula for Cj[k] — e'^ '^°^Cj[k] can be obtained by 
defining 



aj[k] — 2|cj[fc]| cos{ujj[k] + 2^uJoh) 



and 



bj[k] ~ — 2|cj[fc]| sin(a;j[A:] + 2-'cJo^)- 
Now consider the functions /i and /2, given by 



Observe that 



^ |e^2^-o"c,-[fc]-c^^[fc]|2^i||A-/'^j|2 



+ 7II/2-/ 



h\\2 



j,fcez 



/l + /2 _ ,h 



h-h 



where the last equality is a consequence of the parallelogram-law. Since 



j,kei, 



J2 \cj[k]\sm{cjj[k]+V^oh)^P'j^f, 



j.kei 



Y^ \cj[k]\ Ej^k [w CiiOjlk] + 2'ujoh)] 



and 



f' = J2 I9 [fc] |Sj,fe [w{- + Vh) C{lOj [k] + Vuj^h)] 
j,kei. 



we obtain that 



■A + /2 _ rh 



= ei- 



On the other hand, we have the relations 



2 



(./l - /2) = J2 ^''M ^3.k [W C{-UJj\k] - 2^LOoh)\ 

jMez 



and 



lifi - 12) = E l^^[^]l -i^>^ H- + 2''^) C(-^. W - 2-''^o/i)] 



j,kez 



12 



for 



fi - E 4[^]V'..fc and /^' = E ^'l^JV'^.fe. 
Note that fi — f2 — f^ by definition. Hence, we can conclude that 



/1-/2 



/i - h fi' - /a" 



£2- 



This proves the theorem for J = if = Z. 

For the general case, we replace the previous definitions of /i, J2, fi and /j' by 



jeJ.keK 

fi= E «'w^^-.'^ 



/2 = E ^jW'/'i,fc; 



ieJ,/cG-K 






jGJ.fcSK 




similar calculation as before shows that 




E |e'^^-''c,W-c^^[fc]p 




jeJ,keK 




-\\\fl-f^\\l + l\\f2-f^\\l 




1 

~ 2 


h + h fi' + fi' 
2 2 


2 

1 

+ 2 
2 


/l - /2 /i' - 

2 ^ 


3 


_± 


+ ei 









D 



Proposition 14.21 confirms the importance of a smooth localization window w in order to min- 
imize the DT-CWT shift error as previously indicated by Proposition 13.31 

One way to measure the oscillatory behavior of w is by using the concept of Lipschitz 
continuity. By definition, w is called an ^-Lipschitz function if and only if \'w{y) — 'w{x)\ < 
£\y — x\ for every cc,?/ G R. This leads us to the main theorem of this section, which is an 
immediate consequence of Proposition 14.21 



Theorem 4.3 (Decaying rate of the DT-CWT shift error). Let f be a real-valued function 
in i2(M) with DT-C WT coefficients Cj [k] . Consider a translate f^ ^ f{- + h) with DT-C WT 
coefficients c^[fc]. If the localization window w is a compactly supported (-Lipschitz function, 
so that w — w{- -\- 2^ h) is zero outside [p, q], then 



\e''^''^°''cj[k]-c'}[k]\ 



\hcM 



< 2H{q - p) 



(12) 



for every j,k G Z. 



Notice the formal analogy between the DT-CWT bound in (fT2|) and the DWT translation 
sensitivity constraint in (fT|). 
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5 Conclusion 

In this paper, we quantitatively investigated the shift error of the modulated dual-tree complex 
wavelet transform, which can be significantly reduced by performing a phase-compensation 
on the coefficients. 

By introducing a formal description for the DWT translation sensitivity in Section 3, 
we were able to relate the phase-compensated shift error to the shift errors of the real and 
imaginary wavelet components. This study revealed that the superiority of the DT-CWT is 
attributed to the high overall frequency and the slowly varying localization window of the 
DT-CWT. The improved shiftability is particularly noticeable for significant coefficients. 

In Section 4, we estimated the decaying rate of the phase-compensated shift error in case 
of orthonormal and modulated wavelet systems. This allows us to describe the prediction 
error in a similar way as the formal description of the DWT translation sensitivity. 
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